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We study a generalised Gross-Pitaevskii equation describing a d-dimensional harmonic trapped 
(with trap frequency u>q) weakly interacting Bose gas with a non-linearity of order (2 k + 1) and 
scaling exponent (n) of the interaction potential. Using the time-dependent variational analysis, we 
explicitly show that for a particular combination of n, k and d, the generalised GP equation has the 
universal monopole oscillation frequency 2uo- We also find that the time-evolution of the width can 
be described universally by the same Hill's equation if the system satisfy that particular combination. 
We also obtain the condition for the exact self-similar solutions of the Gross-Pitaevskii equation. 
As an application, we discuss low dimensional trapped Bose condensate state and Calogero model. 
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CS| '. I. INTRODUCTION 

> 

There has been renewed interest in the Bose-Einstein condensation (BEC) after its experimental achievement 
by several groups [jlj. The Gross-Pitaevskii (GP) equation explains most of the properties of these trapped dilute 
£S| ■ interacting Bose particles Jl . After the discovery of BEC in a trapped alkali atom, the influences of the dimension 
7— I ' of a Bose system is being the subject of extensive study Recently, low dimensional BEC has been observed in 

MIT Q . It is well known that the nonlinear GP equation for a 2D trapped Bose system has a monopole mode with a 
universal frequency 2cjq HH[rJ ■ ^ n trapped interacting Bose system, the Fermi pseudopotential interaction does 
not affect the monopole mode frequency. This peculiar behavior of monopole mode frequency depends on the scaling 
property of the interaction potential and the order of non-linearity. Without trap potential, the GP equation is scale 
invariant. Introduction of harmonic potential breaks the scale invariance. Pitaevskii and Rosch Q has pointed out 
that due to special property of the harmonic potential, SO(2,l) symmetry still exists. They have shown the universal 
nature of the monopole mode frequency only for 2D GP equation by constructing the generators of the SO(2,l) 
symmetry ||. Here, using the time-dependent variational approach, we obtain a condition for the universality of 
O ■ monopole mode frequency and the dynamics of width of a class of GP equation describing the trapped interacting 
Bose system, at varying spatial dimensionality, order of the nonlinearity and the scaling exponent of the interaction 
potential. Interestingly, the dynamics of the width of a class of GP equation can be described by the same Hill's 
equation if the Bose system satisfy that particular combination or shows the universal monopole mode frequency. 

This paper is organised as follows. In Sec. II, using the time-dependent variational analysis, we show that for a 
particular combination of n, k and d the generalised GP equation describing the <i-dimensional harmonic trapped Bose 
gas exhibits the universal monopole oscillation frequency 2 ujq. We also derive the time-evolution of the width of the 
e?-dimensional trapped Bose system which satisfies that particular combination. The time evolution of the width can 
be described universally by the same Hill's equation which is analytically solvable. We also discuss the condition for 
exact self-similar solutions of the GP equation. In Sec. Ill, we consider quasi-2D trapped Bose gas interacting through 
the Fermi pseudopotential which can be described by an ordinary GP equation. In Sec. IV, we also consider the ID 
modified GP equation which describes the BEC state. We obtain the universal monopole mode frequency and time 
evolution of the condensate width of this system in a time-independent trap as well as a time-dependent trap. In Sec. 
V, we discuss the Calogero type potential. In Sec. VI, we present the summary and conclusions of our work. 
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II. TIME-DEPENDENT VARIATIONAL ANALYSIS 



To show the universality of the monopole mode frequency and the dynamics of the width of a class of GP equation 
(also known as nonlinear Schrodinger equation (NLSE)), we use the time dependent variational method. This vari- 
ational method has been extensively and successfully used to explain many properties in BEC and other nonlinear 
systems. Using this method, Garcia et al. |Tl| derived a general dispersion relation for the low energy excitations of 
a 3D deformed trapped Bose gas which exactly coincides with the results obtained by Stringari fll2| for large number 
of particles. 

Consider a system of Bose particles in d-dimensional harmonic trapped potential Vt(r) = |Mw§r 2 at zero temper- 
ature and interacting with a translationally invariant interaction potential Vj = gdV(r — r"), where r 2 = x\ and 
gd is the coupling constant in d dimension. M is the atomic mass and r is a d-dimensional vector. Under the scale 
transformation, r — ► Ar, we assume Vj(Ar) = X n Vi(f) where n is the scaling exponent of the interaction potential. 

We are considering a generalized d-dimensional GP equation E3 | describing a d -dimensional harmonic trapped 
weakly interacting Bose gas with a non-linearity of order (2k + 1) and scaling exponent n of an arbitrary interaction 
potential Vj = gdV(r — r 1 ) 

*^7^ = l-^ 2 + V t (r) + g d [ d m r'r(r )V(r- ^(^(r)^ 1 ^]^, 



ih ot = h 2M v + Vt{r) +9d J ^WMr- fMnrnvr* WW, (!) 

where V 2 is the d-dimensional Laplacian operator and * denotes complex conjugation. The generalised d-dimensional 
energy functional is 



(2) 



If m = d = arbitrary dimension and V(r — r 1 ) is a d-dimensional delta function potential, the above equation gives 
the non- linear term of order (2k + 2). If k = 1 and m = d > 1, Eq. (Q) becomes ordinary GP equation, and if k = 2 
and m = d = 1, Eq. (|l|) gives the ID modified GP equation in which the \ip\ 6 interaction term is present [pT| , p^[ . 

One can write down the generalised Lagrangian density corresponding to this system which is 

L = Y^^T ~ r ^ ) + [ ^ V ^ 2 + ^l^l 2 + (fcTT) / dm r'r(r)r(r")V(r- ^)^mr)\^^(r)]. 

(3) 

Without interaction potential, Eq. (0) reduces to the linear Schrodinger equation which has the Gaussian wave 
function for the ground state. In order to get the evolution of the square of the radius, we assume the following 
Gaussian wave function, 

i>(r) = C(t)e 2a « a(t) , (4) 

where C(t) is the normalization constant which can be determined by the normalization condition J \ip\ 2 d d r = N 
and cio = yJh/MuQ is the oscillator length. N is the total number of particles in the system. a(t) and (3(t) are the 
time-dependent dimensionless variational parameters. a(t) is the square of the radius of the system. In order to 
describe the dynamics of the variational parameter a(t), the phase factor i(3(t)r 2 is required Note that the same 
kind of phase factor is also present in the Eq. (4) of the Ref . J 1Q|] to get the dynamics of the width of a system. We 
get the effective Lagrangian L e ff(a,/3) by substituting Eq. ([I) into Eq. and integrating the Lagrangian density 
over the space coordinates, 

Leff = ^H[- a 0d + - + ad(l + (3 2 ) + G d a^ +m - kd ^ 2 }, (5) 
4 a 
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where 

G d = G J d d Rd m R'e-( R2+R ' 2 )V(R-R')e-( R ' 2 + R2 ( 2k - 1 », (6) 
where Gn - 4 g^ 2(d+T " + " )/ M" + ""'"" and r 2 _ 2aR 2 a 2 

Using the Lagrange equation of motion, we get the equation of motion for the variational parameter a(t). The time 
evolution of the variational parameter a(t) is 

& = 2a/? 2 + - - 2a - {n + ™ ~ k ^ Gd a (n+m-kd)/2 (7) 
a a 

The " ." represents the derivative with respect to the dimensionless variable r = uj^t. The total energy of the system 
is 

E = 4E/Nhu = af3 2 d +-+ad + G d a {n+m - kd)/2 . (8) 

a 

At equilibrium, (3 = and ^| a=ao = 0. The equilibrium point «o satisfies the following relation in c£-dimension 

2 (n + m- kd)G d i n +m-kd)+i 
a + — a' = 1. (9) 

We choose n, m, k and d such that they satisfy the following relation 

n + m- kd +2 = 0. (10) 

If n, m, k, and d satisfy the above relation, Eq. (Q) can be written as 

a + 4a = 2E(a). (11) 

Define I(t) = J d d r\ip(r, t)\ 2 r 2 which is the expectation value of the square of the radius. The above equation can 
be re-written as 

From now onwards " ." indicates derivative with respect to the time t. For monopole excitation, the average value of 
the collective coordinate r 2 oscillates around its equilibrium value. We expand I(t) around its equilibrium point Iq 
as I{t) = Jo + SI(t) and we get, 

5I + 4u 2 SI=^, (13) 

where Eq = d{-^- + ao + ^f) is the equilibrium energy of the system. This equation is identical with the classical 
equation (1) of Ref. Q which was derived using the virial theorem. The solution of this equation is 



I(t) =Acos(2u> t + 6) + — V (14) 

Mbjf 



o 



A and are constants which can be determined from the initial conditions on ?/>. Thus the square of the system 
radius oscillates with a frequency 2wo in any dimension if n + m — kd + 2 = 0. This monopole frequency 2 loq is 
universal because it is independent of the interatomic interaction, dimensionality of the trap, order of non-linearity 
and the total number of atoms N. For ordinary GP equation in c?(= m) dimension, k = 1. Hence the monopole 
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mode frequency of an ordinary d-dimensional GP equation is the universal frequency 2u>o if n = —2. Eqs. (0) and 
( |To| ) explicitly shows how the scaling exponent of the interaction potential and the order of non-linearity gives the 
universal monopole mode frequency. 

If n, m, d and k satisfy the Eq. ([To|), the equation (0) can be re-written for time-independent trap as well as for 
time-dependent trap as 

^~m + ^ {t)m = m (15) 

where Q d = (^) 2 (1 + G d ) is the positive invariants under time evolution. This Q d does not depend on the frequency 
of the trap potential. Let X(t) = y/ I(t) is the condensate width. The above equation can be written as 

X^+^Xit)^^.. (16) 

This is a non-linear singular Hill's equation in d dimension. This is an interesting result because the time evolution 
of the width can be described by the same Hill's equation in any dimension ifn + m — fed + 2 = 0. The dynamics of 
the width is universal because the time evolution of the width can be described universally by the same non-linear 
singular Hill's equation in any dimension. Moreover, this Hill's equation can be solved analytically and exactly. The 
Eq. ( |l6| ) can be viewed as the classical motion of a fictitious particle in an effective d-dimensional potential 

_ u>j(t)X 2 Q d 
eff ~ 2 2X 2 ' 



The general solution of Eq.(|l6|) is X(t) = \Ju 2 {t) + ^jfzv 2 {t) where u(t) and v(t) are two linearly independent 

solutions of the equation p + ui 2 (t)p — which satisfy the initial conditions u(t ) — X(to), u(to) = X'(to), v(to) = 
and v'(to) ^ 0. W is the Wronskian. This time-dependent Hill's equation can be solved explicitly only for particular 
choices of u>o(t). For time-dependent periodic trap potential, the equation is known as Mathieu equation which is 
well studied. In most of the BEC experiment, the trap potential is time- independent. For time-independent trap 
potential, a solution of this equation is 



X(t) = Jcos 2 (2uj a t) + -^4-sin 2 (2w t). (18) 



From Eq.(|To|), we get kd = 2 for ID modified GP equation, and 2D ordinary GP equation with n = —2. This is the 
condition for self-similar solution of a non-linear GP equation discussed in by Rybin |l4|,[l5| . There is no self-similar 
solution in ID ordinary GP equation since it does not obey Eq.(pX|). 

We now discuss few systems which satisfy the Eq. ( fTo| ) and shows the universal monopole mode frequency and the 
dynamics of the width of the system that can be described by the Eq. (J18|) . 



III. 2D TRAPPED BOSE GAS 

It has been shown by Baganato et al. [|l6| that for an ideal 2D Bose gas under harmonic trap a macroscopic 
occupation of the ground state can exist at temperature T < T c = ^N/((2)hu/kB- With present technology one 
can freeze the motion of the trapped particles in one direction to create a quasi-2D Bose gas. In the frozen direction 
the particles execute the zero point motion. To achieve this quasi-2D system, the frequency in the frozen direction 
should be much larger than the frequency in the X-Y plane and the mean interactions between the particles. The 
wave function in the z-direction is 
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ip(z) = Ce (19) 

where a z = yjfi/muj z is the oscillator length in the z-direction and C is the normalization constant. When we integrate 
out the z-component in the 3D GP energy functional, then we get the effective interaction potential Vi = g2^ 2 (r — f ) 
where gi = 2\/2~nfiuj z aa z is the effective coupling strength in 2D. a is the s-wave scattering length in 3D. This 
quasi-2D BEC state can be described by the ordinary GP equation which is valid only when a < a z j5|. The same 
effective coupling constant is obtained in Under the scale transformation, r — > Xr, Vj{Xr) = j?Vi(r). So this 
interaction potential has the scaling exponent n — — 2. We will show that for an isotropic harmonic trapped system, 
the monopole mode frequency is 2dJo- In this system, m = d = 2 and k = 1 which satisfy the Eq. (|l0|). 
The time-evolution of this system is 

X + ^pf = !f. (20) 

The same equation is obtained by using different method at zero temperature |jl^,0. Here, we have identified the 
value of Qi at T = which is Q2 — N M 1 (1 + yl^N-^-). The dynamics of the width can be described by the Eq. 
(H). At equilibrium, X$ = %. Whe n we linearize Eq. (^o|) around the equilibrium point Xq, we get, 



5X + Alu^SX = 0. (21) 

Once again we obtain the universal oscillation frequency of the condensate width to be w = 2loq, corresponding to the 
frequency of a single particle excitation in the condensate state. 



IV. ID TRAPPED BOSE GAS 

It has been shown by Kolomeisky et al. jL5| that the GP equation describing BEC in ID contains |i/'| 6 interaction 
term instead of familiar \ip\ 4 term. The energy functional of ID dilute interacting Bose gas at zero temperature is 

In this system, d = m = 1, k = 2 which satisfy the Eq. ( |l0| ) and g\ — ^jj-- 
The time-evolution of the system is 



*+ w o*=3p- (23) 



where Q\ = ^1 [1 + ]■ The same equation is obtained in [J14 15 by using different method. But here we 
have identified the value of constant Q\. This ID modified GP equation also exhibits an universal monopole mode 
frequency 2tt>o and satisfy the self-similar condition. The dynamics can be described by the Hill's equation which is 
exactly solvable. 



For ID ordinary GP equation, m = d = k = 1, which does not satisfy the relation (10 ). So the monopole mode 
frequency is affected by the interatomic strength and number of particles and the time-evolution of the system can 
not be described by non-linear singular Hill's equation. 
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V. CALOGERO MODEL 



Only an inverse square pair potential, 1/r 2 , has the scaling, n = —2, in any dimension and it satisfies the Eq. (|l(]| ) 
for k — 1 in arbitrary d. A system of trapped particles interacting through pair potential 1 /r 2 is well known as the 
Calogero-Sutherland model fig] . In Ref. |19j, they have shown that the energy spectrum of a system of harmonic 
trapped particles interacting with 1/r 2 potential is divided into sets of equidistant levels with separation 2ojq in any 
dimension. In Ref. fl2(i|| , they have considered a one-dimensional system of particles interacting by an inverse square 
pair potential and shown that the monopole mode frequency is the universal frequency 2luq. For m = k = and 
n = — 2, we get Calogero model |Is| in arbitrary d which gives the same universal monopole mode frequency 2loq. In 
this system, the dynamics of the width can be described universally by the Eq. (|l§|). 

Apart from these two potentials, there are some other model potentials which have the scaling exponent n = —2. 
For example, the derivative delta function potential in ID, the double delta function potential in ID, and the 1/x 
delta function potential in ID. For all these model potentials, one can get the universal monopole mode frequency 
and the dynamics of the width. 



VI. SUMMARY AND CONCLUSIONS 



In this work, using the time-dependent variational analysis, we have obtained the condition, n + m — kd + 2 = 0, 
for which a class of NLSE exhibits the universal monopole oscillation frequency 2 lj . This monopole mode frequency 
is universal because it does not depend on the strength of the interaction potential, dimensionality of the trap, order 
of non-linearity and total number of atoms N. We have also shown the universality of time evolution of the width in a 
time-independent trap as well as a time-dependent trap for a class of GP equation. Interestingly, this dynamics of the 
width can be described by the same Hill's equation ( |l6| ) if the Bose system satisfy the condition, n + m — kd + 2 = 0. 
This Hill's equation is analytically and exactly solvable. Using the time dependent variational analysis, we have also 
obtained the condition for self-similar solutions, kd = 2, of a class of non-linear GP equation. 

As an example, we have shown that the 2D trapped Bose gas describing by the ordinary GP equation, ID BEC 
describing by the modified GP equation and a system describing by the Calogero model shows the universal nature 
of the monopole mode frequency and the dynamics of the width of the system. We have identified the values of the 
positive invariants quantity, Q2 and Q\ f or 2D and ID BEC. 

After submitting this paper on the net ( |cond-mat/001218S ) , there is another paper H] appeared where the univer- 
sality of the monopole mode and the dynamics of the width of a class of NLSE has been shown exactly and analytically 
which is obtained from different method. Note that the phase factor in the variational ansatz of the wave function 
(eq.(4)) is required to get the dynamics of the widt h [O . Similarly, the phase factor is also present in equation of the 
time dependent transformations in 1 10 2^]. In Ref. ]2l||, using those time dependent transformations, the universality 
of the monopole mode and the dynamics of the width is obtained for a class of NLSE. Our results obtained from the 
variational approach are identical to the exact results of pl[ | . This is because of the similarity between the variational 
ansatz for the wave function and the transformation equation (eq.(2)) in Ref. 21 1. 

Our result is also valid for SO(2,l) invariant multicomponent NLSE describing BEC state. Recently, 2D and ID 
BEC has been observed in MIT Q . Experimentally, it is also possible to check the validity of the two-body potential 
in 2D BEC state || and the order of the nonlinearity in ID BEC state by measuring the universal monopole mode 
frequency and the dynamics of the width of the condensate state. 

I would like to thank Pijush K. Ghosh, M. V. N. Murthy and K. Santosh for helpful discussions and many valuable 
suggestions. 
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